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SUMMARY 


A simplified method has heen developed for the calculation of heat 
transfer and skin friction in the compressible laminar boundary layer 
with an arbitrary Prandtl number near unity and an arhitiary strearawlse 
pressure gradient and wall tea^jerature distribution. By comparison of 
numerical resvilts with some exact solutions, the method is shown to 
give accurate results for the case of houndary-layer cooling when a 
certain fifth-degree polynomial is used for the thermal profile. The 
use of this polynomial also results in recovery factors which are 
accurate to within a few percent on both the flat plate and the circu- 
lar cylinder normal to the flow. 

The method may he extended to calculation of heat transfer with 
equilibrium dissociation provided the wall temperature' is below dis- 
sociation ten^eratiires . A nvunerical exan^jle for the heat transfer in 
the stagnation region of a blunt body shows that equilibrium dissocia- 
tion probably has little effect on the actual heat-transfer rate. 

The effect of lateral radius of curvature on skin friction and 
heat transfer is investigated for con 5 )ressible flow. The approximate 
method -used herein indicates that this effect is negligible unless the 
square root of the 'wall Reynolds number (based on body length and the 
gas properties e'valuated at the waU) is of the order of or less than 
100 times the body fineness ratio. 


INTRODUCTION 


The calculation of the laminar boundary layer -with an arbitrary 
pressure gradient is as yet ln 5 )ractlcal except by the Von Karman- 
Poblhaiisen integral methods . In general, these methods are somewhat 
tedious to apply; however, recent in 5 )rovements have reduced the calcu- 
lations -with zero heat transfer to a simple quadrature of certain 
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functions of the known free-stream velocity and temperature distribu- 
tions. These various refinements which both sin^jlify and improve the 
accuracy of the basic integral methods have been developed by Thwaites 
(ref. 1), Truckenbrodt (ref. 2), and others. (See discussion in ref. 3*) 
Rott and Crabtree (ref. 3) have extended the method of Thwaites to com- 
pressible flow with zero heat transfer and a Prandtl nuniber of unity. 
Morduchow and Clarke (ref. 4) have also treated the problem of con5)res- 
slble flow with zero heat transfer and a Prandtl number of tinity by a 
different approach which, for practical purposes, gives results of equal 
accxiracy aufl equal sin5»licity of calculation. The method of reference 4 
results in velocity profiles of sufficient accuracy for stability calcu- 
lations, whereas the methods of references 1 and 3 cannot be used for 
this purpose. 

If the heat transfer is included as another unknown, the problem 
is more complicated, particularly for an arbitrary surface ten5)erature 
distribution. Dienemann (ref. 5) shown that the calculation with 
heat transfer in Incompressible flow can be sin5)lified by application 
of the Holsteln-Bohlen version of the Von Karman-Pohlhausen method. 

For an Isothennal surface and a Pi'andtl number of unity, Morduchow 
(ref. 6) uses assiui5)tions similar to those in reference 4 to simplify 
the problem of heat transfer on a sweat-cooled wall in con5)ressible 
flow. The work of reference 6 is an extension (to include a normal 
velocity at the wall) of the method of Kalikhman (ref. 7) which also 
gives detailed results only for a Prandtl nximber of \inity but is not 
restricted to an Isothermal surface. The stagnation-temperatirre pro- 
files are assumed as fourth-degree polynomials in the normal distance 
in both references 6 and 7 • 

Recent theoretical and experimental work (for example, refs. 8 to 
11) bRs indicated that boundary layers in high-speed flow may be la m i n ar 
at relatively large Reynolds nymbers under conditions of boundary-layer 
cooling and favorable pressvire greidients. Therefore, the accurate cal- 
culation of heat transfer and skin friction under these conditions is 
of considerable Interest. The purpose of the present paper is to trace 
the development of a method by which these calcTilations can be made 
with relative slD5)llclty. The present method, which is essentially a 
modification of Kalikhman 's method (ref. 7), uses an arbitrary Prandtl 
number near unity and several different stagnation-tenqoerature profiles 
in conjunction with a fourth-degree velocity profile. The principal 
siD5)lifylng assun^itlons and procedures are (l) equal thermal and veloc- 
ity boundary-layer thicknesses, (2) linear viscosity- teu^jerature rela- 
tion, (3) use of the first coefficient in the polynomial for the 
stagnation-ten^jerature profile as one of the unknowns in the final solu- 
tion of the momentum and energy equations, and (II-) application of the 
Holstein-Bohlen method to avoid the use of the second derivative of the 
stream velocity. The results in the form of heat -transfer an d skin- 
friction coefficients anfl recovery factors are then compared with exact 
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solutions for the nonl so thermal flat plate (ref. 12 ), stagnation flow 
(refs. 15 and 1 ^), and flow over a cylinder (ref. 15 )- The recovery 
factors on a cylinder are compared with the supersonic data of Eher 
(ref. 16). 

The effects of equilibrium dissociation (as given by Moore for the 
flat plate in ref . I7 ) on the heat transfer at a surface cooler than 
about 2,000° H in the presence of a velocity gradient are also considered 
in the present analysis . A numerical example is then given that shows 
the approximate effect of dissociation in the stagnation region of a 
blunt -nose body. 

After the present investigation was con5)leted, the author learned 
of a paper by Paul A. Libby and Morris Morduchow of the Polytechnic 
Institute of Brooklyn who, in some instances, vised procedures and 
assvm^tlons similar to those of the present analysis. Their method 
used a sixth-degree polynomial for the velocity profile and a seventh- 
degree polynomial for the stagnation-tengierature profile. Althou^ 
the sixth-degree velocity profile is probably necessary for stability 
calculations, the seventh-degree thennal profile, as indicated by this 
analysis, in some cases gives less accurate results than a singiler 
fourth- or fifth-degree thermal profile. 


SYMBOLS 


Ai,A2, 

tii,h2, 

Cf 


coefficients of velocity profile (eq. (l 2 )) 
ajj constants in expression for wall teng)erature (eq. (47)) 

coefficients for stagnation-enthalpy profile t (eq. (16)) 
coefficients for the density function (eq. (26)) 
mean skin- friction coefficient 


Cw - 


^VPw 

tilPl 


Cf 


°P 

Fr 


local skin-friction coefficient 
specific heat at constant pressure 
"local" recovery factor (eq. (62)) 

"local free stream" recovery factor (eq. (63)) 
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KojK2_, . • • 

k 
L 
M 


m 


Wu 

n 


P 

<1 





^ PqUqL 

Po 

^ PlU^L 

Pi 



Poo^J- 


stagnation enthalpy. 



static enthalpy 
constants (eq. (75)) 

conductivity 
characteristic dimension 
Mach number 

constant, 

Niisselt nvpnher 

degree of polynomial for wall ten^jerature (eq. (il-y)) 
pressure 

local convective heat-transfer rate 


r 


T 

Te 


radial distance from axis for cylindrical coordinates 
absolute temperature 

wall ten5)erature for zero heat transfer 
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T = 

^T.r 


T, 


w 






t 

* 


stagnation temperature, T + 


u‘ 


2Cr 


t* H - % 

t = 


L* 


2_ % ' 


t* = H - Hw = H - Cp^Tw for < 2,000 R 

u velocity in x-direction 

V velocity in y-direction 

V velocity ratio, u/uj^ 

x,y boundary-layer coordinates 5 y = 0 is body surface 


Yn'(o) 

Z 

OC* y CX^ j ^^^*2 


heat-transfer parameter for nonisothermal flat plate 
in notation of ref. 12 

confutation parameter defined by equation (67 ) 

constants in energy-loss thickness (eq.. (55)); numerical 
values given in table I 


= cos 


-1 ^ ~ 


7 

A 


ratio of specific heats (approximately 1.40 for air) 
nominal boundary-layer thickness in x,t) plane 


A* = A f (1 - t)dTi* (eq. (25)) 
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T 1 


* _ 


a 

A 


displacement thickness in x,ti plane (eq.. (25)) 

ry 

Dorodnitsyn vajriable, / dy 

Jq P 1 

modified Dorodnitsyn variable for slender body (eq.. (AT)) 



0 

0 

A 

A 


0 

L 




moment\mi-loss thickness in x,q plane (eq. ( 22 )) 


modified Pohlhausen parameter (eq. (15)) 

modified Pohlhausen parameter for slender body (eq. (A 14 )) 

viscosity 




, ^ ^ / Pi 

ypi^i^ 


P 



a 

T 

tP 



mass density 
stagnation density 
Prandtl number, Cpp/k 
local sheax stress 

energy-loss thickness integral in x,t) plane (eq. ( 24 )) 
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a 

Subscripts : 
1 
00 

V 

B or o 
P 


density integral. 




local free stream just "outside" boundary layer 
free stream at infinity 
wall values 

some constant reference quantity 
flat-plate values (zero pressure gradient) 


BASIC EQUATIONS 


In the following discussion, the boundary-layer equations and 
their integrals are given in their most general form where the only 
restrictions are that the gas is a continuum and that the boundary 
layer is thin compared with a characteristic dimension of the body, 
that is, the characteristic Reynolds number of the body is large. For 
application to dissociation, two additional restrictions are used. 

These restrictions axe as follows: the assumption of equilibrium dis- 

sociation throughout so that no rate process terms are required in 
the energy equation (see ref. 17 and the restriction that all wall 
tenqjeratures are less than about 2,000^ R. The two-dimensional equations 
only are considered herein since the results may be extended to axlsym- 
metrlc flow by Mangier 's transformation (ref. I 8 ) except where conditions 
are such that the bovmdaiy- layer thickness is of the order of the lateral 
radius of the body. The effects of these conditions are investigated in 
detail in the appendix for a slender cylinder alined parallel to the 
flow. 


The boundary- layer equations for steady flow are written as follows: 


Momentum equation: 


pu 




^ 8 / 8u\ 


(1) 


0 


_ ^ 


( 2 ) 
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Continiiity equation: 


J- pu + ^ py = 0 


( 3 ) 


Energy equation: 


8 h 

pu ^ + pv — = u 


Sy \ V3y/ 


ih) 


or in terms of the stagnation enthalpy H = h + — , 


s * i ' |r(‘‘ i) ^ 


A 

P - 1^ 

/ ^ 

By\0 / 


( 5 ) 


where the Prandtl nimiber 0 is defined as 




- U <3^1 

k dr 


Integration of equation (l) across the boundary layer (with the 
conditions u = v = 0 at y = 0 and u = U]_ at y = ~) and the use 
of equations (2) and (5) result in 



dn ‘ 3 -Ut 

where Bernoulli's equation 3^ = “Pl^l j — 3 ias been used for the pres- 


dx 

sure gradient, and the free -stream stagnation density 
constant along the edge of the boundary layer. 


p*^ is assumed 
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Similarly, integration of equation (5) and the use of equation (3) 
result in 


_d^ 
dx , 


'0 


-1 - 




dy + 


1 ^^1 


1 r u 


^^1 dx t*3^ ^ 0 


'1 - 


h*, n* 


dy = 


(h t* \ 


(T) 


where the free-streain stagnation enthalpy is assumed constant. 

Transformation to the Dorodnitssm variable tj defined as 


T] = 




( 8 ) 


results in 



This derivation of equations (9) and (lO) follows that of reference 7. 
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VELOCm AHD EHTHALPy PROFUiES 


The usiial procedure in the Von Kaniian-Pohlhausen integral methods 
is to assiime the velocity profile as a polynomial in the normal distance 
parameter t]. The fourth-degree polynomial is sufficiently accurate for 
most engineering purposes except stability investigations and prediction 
of sepaxatlon as shown, for exan^ile, in references 4 and 19 for zero 
heat transfer . There is no way of knowing beforehand whether the same 
conclusions are valid in the presence of heat transfer- However, since 
the present analysis is concerned mainly with the effect on the heat 
transfer of different enthalpy or teii5)erature profiles, the fourth- 
degree polynomial for the velocity profile is used throughout for 
6in5)licity in calculation. 

The velocity profile is then written 

^ = w = An* + A2 (t}*)^ + + A4(t]*)^ (11) 

The usual hoimdary values are specified at the edge of the boundary 
layer such that, at -q* = 1 , 


w = 1 


bv 




= 0 


From these boundary conditions and equation (u), the following relations 
are obtained; 


Al 


bvr \ ^ 12 + A 


^3 - I 




3 (q*)- 


5A - 12 






(12) 





A = -2A2 = - 
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In order to sln^llfy the remaining houndary condition, it is convenient 
to introduce the asstm^jtions that, in the immediate vicinity of the 
surface only. 


K = ? T (13) 

■‘■w 

and 

P = Pw^v ^ (1^) 


These assun^jtions should not he unduly restrictive so long as the surface 
ten^ieratiire is less than 2 , 000 ° R (that is, well helow the dissociation 
temperature of the gas) since may he evaluated as accurately as 

desired to correspond to Ty. Evaluating equation (l) at y = 0 n.nfl 

using equations (13)3 (l^)j ( 8 ), and (l 2 ) give the required expression 
for A as 




Ha" 

dx 


(15) 


In a similar manner the stagnation enthalpy profile is written as 
a general polynomial in q*. 


^ = t = Bpq* + + B)^(q*)^ + . . . (16) 

where up to six terms will he used. The edge boundary conditions are 
taken as two or more asyn^itotic requirements at q* = 1 and are 


t = 



8(q*)^ 


0 

J 


(IT) 


where n is a positive integer. 
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Equations (17) in5>ly the assun^tion that the thermal and velocity boundary 
layers are always' of equal thickness. In contrast to this assun^ption, 
Morduchow (ref. 6) has assvuned that the thermal and velocity layers are 
always in the same ratio on any particular body. The results of the 
exact solution of Bro'wn and Donoughe (ref. I5) of wedge-type flows 'vrf.th 
Prandtl number of O.7, however, indicate that the nominal value of this 
ratio is. approximately uni'ty for a ■wide range of pressvire gradients and 
■wall temperatures. Fiir thermore , the computations of Le'vyN^ (ref . 20 ) 
show that, e-ven for a Prandtl number as large as the thermal boundary- 
layer thickness can be considered an roughly JO percent of ^the ■velocity 
boundaiy- layer thickness because of the asyn^jtotlc nature of the profiles. 

As ■will be sho'wn subsequently, the resiilts of the present analysis indi- 
cate that the assun^itlon of equations (17) gives good accuracy when used 
in a cer-bain combination ■with the wall botindary conditions for equation (16). 


These ■wall boundary conditions are deteimlned from eqimtlon (5) 
and its first deri-vative -trLth respect to y e'val u ated at y = 0 . With I 
the use of equations (15) and (l 4 ) and the transformation - formulas fr'om 
equation (8), the ■wall boundary conditions for the qhthalpy profile may 
be ■written as . ■ 



where the additional assun5)tions ha've been used that o and Cp are 
constant very near the ■wall. By using equations (l 2 ) and the assuii5)tlon 
that the thermal and ■veloci'ty boundary layers are of equal thlckuess, 
equations (18) and (19) can be written as 


2B2 


8^ 




(20) 



8 ^ 

8(ti*)^ 


A 

t\ 6 




^ dTy/dx 

pp duq^dx 


3(1 - (%)uq 


(Eq. (15) has also been vised in eq. (l 9)0 


( 21 ) 
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GEaiEHAL METHOD OF SOLUTION 


The general procedure for the solution of equations (9) and (lO) 
consists of using two or more of the asyngitotic requirements (eqs. (17)) 
and one or hoth of the wall conditions (eqs. (20) and (2l)) such as to 
specify all hut one of the coefficients of the stagnation enthalpy pro- 
file in terms of A and the known external flow. The remaining coef- 
ficient of equation (16), say Bq, and A will then he the two unicnown 
quantities to he found from the simultaneous solution of equations (9) 
and ( 10 ) . 

Before the solution can he effected, explicit expressions for the 
integrals of the velocity and enthalpy profiles are required. Thus, 
as usual for the fourth-degree velocity profile, the momenttim-loss 
thickness integral is 

I = (w - . Uj(37 - I - 2g) (22) 


and the displacement thickness is 



36 - A 
120 


(25) 


The general expression for the energy-loss thickness integral is 
obtained as 


9 

A 


w(l - t)dTi* 


10 30 kp ^ pfin 0 ^ P ppn ° 


1+2 


280 


36 


630 5 " 220 



Bp - ^ 


(2l+) 
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where ^ 3 ^ • • • will he expressed in terms of throiogh the 
houndary conditions (eqs. ( 1 T)j ( 20 ), and ( 2 l)). The corresponding 
expression for A*/A is 


A 



% 

5 


6 


7 


(25) 


which will he used in the solution for the flow of a perfect gas. 

The only integral in equations (9) aiid (lO) as yet iinspecified 
contains the density profile Pi/p. For the general case where equl- 
lihrium dissociation will he considered, the density ratio as determined 
from the data of reference.17 can he sp,tisfactorlly approximated hy a 
fifth-degree polynomial in the enthalpy h and is 


£1 = £1 

P Pb 



^2 4 + ^3 






(26) 


In general, the quantities hQ_, h2, . . . will depend on the local 
pressure; however, as a first approximation these quantities are treated 
as constants. For a perfect gas and constant Ct>, equation (26) reduces 
to 


= A = JL 

p hi Ti 


The enthalpy ratio h/hg is expressed in terms of the stagnatlon-enthalpy- 
drop profile ^ hy the relation 

A = 

hB hg 2hg \^hg hB 2hB/ 

from the definition of t. For a perfect gas, this equation reduces to 


1^1 


Ti 2hi 





MCA TK 5005 


15 


The density integral 




can now he expressed 


in terms of A and the stagnation enthalpy profile coefficients 
Bp, .... This integral is too lengthy to repeat here for the 
general case with dissociation; however, for a perfect gas, it becomes 


ft 

A 




(27) 


The local skin-friction coefficient is by definition 



or, when transfonned to t]*. 


Cf = 


^^w ^w i/^vr \ 


(28) 


Then, by using ecpiations (I5) and ( 12 ), the general relation is obtained 


Cf 


^ ./Py X ^ui~ 12 + A 
1 /Pl ui dx A 6 


(29) 


which is dependent only on the wall temperature, the external flow, and 
the parameter A. 

Similarly, the local convective heat-transfer rate is defined as 



or, when transformed to t]*. 


Pw 


Pw / 5 t \ 1 


(50) 
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since | = c-n,/— ] . Use of equations (15) and (16) in equation (50) 

/w 

then gives the relation 


V u-j 


V Cp^VVPw ^ 


(31) 


"which is dependent on the wall tenipfira'turs ^ the exte rn a l flow, and 
hoth 7 \ and 

APPLICATION TO FLOW OF A PERFECT GAS WITH CONSTANT Cp AND a 


General eqiiations .- When the results of the preceding section are 
used for a perfect gas together "with the assumption that the thermal 
and "veloci"ty houndary layers are of equal thickness, equations (9) and 
(10) can he "written as 



As an illustration of the detailed procedure, a fifth-degree tempera- 
ture pniyn mni al is used in the following de"velopment . This polynomial is 
chosen so as to satisfy "two conditions at the edge of the boundary layer 
ant^ two conditions at the "wall. At the edge of the boundary layer, that 
is, at q* = 1 , 


t 





( 54 ) 



RACA TR 5005 


17 


These two conditions at t)* = 1 give 

Bij. = 5 ~ ” 5 B 2 “ 2B^ 

and 

= -^ + 3 Bq_ + 2B2 + Bj 

Equations ( 24 ) and (25) may then he written as 


— = q + + a2B2 + “■3^3 + °6®2 “7^5) ^ 55 ) 


and 


A 3 5 15 60 5 


(56) 


where the a's are constants and B2 and B3 depend on the known 
external flow, the parameter A, and the wall temperature distrihution 
as determined from the two wall conditions (eqs. ( 20 ) and ( 2 l)). Rumerl- 
cal values of the a constants are given in table I. For a perfect gas 
and constant Cp, these mil conditions can be written as 


2L-1 fl2 + A' 


B2 = (1 - cr)- 


Mi' 


T^l/Ti 


i^y 


1 - T 


( 57 ) 


w 


and 




_ a 12 + A a2 (p\) ^1 1 


dT, 


7-1 


w 


PwV^w 


1 _ fp dx 
■‘•w 


- (1 - ct) 


Mt 2 A 


12 + A 


(38) 


T*l/Ti 1 - Tw 


T 


where T^ = -^ 
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du-i du-i 

Conipeirison with flat-plate flow: — ^ = 0.- For — ^ = 0^ the 

^ ^ 

momentum equation (32) reduces to ~ 

^ = 2 

” ui(p\)^A 

since A = 0. By using equation (22) (with A = O)^ this equation can 
he Witten as 


37 dA _ 

315 dx ui(p\)^A 


(59) 


which by integration gives 


2 = 1, 3^ 1 r 


PwPw 


(40) 


Equations (3?) and (38) can now he witten as 


and 


7-1 


B2 = 4(1 - a)- 


Ml^ 


<1 - Tv) 


(41) 


B = 420 „ «Try/dx 1 


57 1 - Tv ^0 


r 


^HrPu- ^ 


(42) 


The energy equation (eq. (33)) is written for the flat plate as 


^ 1 - Tv ou]^(p*-. ) A 


( 43 ) 
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or by using equation (35) -with A = 0 and equations (4l), (42), (39), 
and (40), equation (43) may be written as 

^ + B ^ CT ^ f ^ 

dx ^ 2f^V 630 aaqy 1 - 37 “"I ^ 1 - T^- 


^ . 2^ CT ^(1 ^ 2 . 1 L 4o^(l - CT). 

0-1 57 “I \ ^ - Tw ^ 


Mq^ 

_2 ^ 1 

T*i/Tq 1 - T 


where = — — / jxwPw Integration of equation (44) gives the 

^■wPw Jq 

required solution for Bq* 

In order to conq)are the results directly with those of reference 12 
it is convenient to set f^ = x, whereupon the expression for Bq is 


T, _ 1 (X)-”1 

Jdt — — 

^ Tv - 1' 


% 2 ^ g . ^ ^ - 1 + 

X 7 -lpTv-l 2 

1 + Mq 2 


o^cr - a)::: + — a^ox XT^ - l)5Pdx (45, 

\57 ^ ^Tv - 1 57 ^ Tv - 1 

where x = ^ and m = — - . From equations ( 30 ) and (4o), the 

L 2 126O0U a 


heat transfer is 


= - 1 V) a - ^ ^ 
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where Cy = — ^ = tlaiEit for a constant and cr. The Chapman- 

Ruhesin polynomial for the wall temperature distribution (ref. 12) is 


^w 

Ti 






(i^7) 


Substituting equation (45) (with T^ given by equation (47)) into 
equation (46) and integrating result in 


<lw 


- 1 /37 

^V315 


^1^1^ 





liiO ^2 
57 


+ 




a 

2 


n + m 


(^8) 


Thus, in order to evaluate the accuracy of the heat transfer obtained 
from the present method by using a fifth-degree temperature profile, 
the expression 


2 ^ 

/37 


420 ^ 

/ 210 

\ (X 

ap 

/315 


n*^ + 

37 

(“5" 37 • 

- ajn - - 


n + m 


(49) 


is to be con 5 )ared with the quantiiy -Yjj'(o) in the notation of refer- 
ence 12.. This comparison can be obtained from figure 1 where expres- 
sion ( 49 ) is plotted against n for several different temperature pro- 
files. (Numerical values of the a constants are given in table I for 
several of the thermal profiles used. ) The results are Identified with 
the corresponding temperature profile by designating the degree of the 
polynomleil and the niimber of wall and, edge boxmdary conditions used to 
form the temperatiare profile. The quantity -Yj2’(0) from reference 12 

is plotted in the same figure. Also shown in this figure is a table of 
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the percent error in the heat transfer as computed from the various 

dT 

thermal profiles for n = 0 which corresponds to the case with — — = 0. 

dx 


Con^iarison of the resiilts of the present method with the Chapman- 
Ruhesin nonlsothermal flat-plate solution shows that the sixth-degree 
(1 wall boundary condition, k edge boundary conditions), fifth-degree 
(2 wall, 2 edge), and sixth-degree (2 wall, 3 edge) stagnation tempera- 
ture profiles may be expected to yield satisfactory heat-transfer com- 
putations when the wa ll temperature distribution can be approximated 
by a polynomial in which the lower order terms (n < 5) predominate. 
However, it is in^jortant to note that a polynomial for the wall tempera- 
ture is not required in the present method. Thus, a general conclusion 
applicable here is that the heat-transfer values would be in some error 

towards the rear of the body ^x > ^ if the temperature distribution 

in this region is characterized by large gradients . The reason for 
the error in the present method under this condition is probably 
the assuD^jtion of equal thermal and velocity boundary- layer thicknesses 
since the exact solution of reference 20 shows that the thickness of 
the thermal layer is reduced by increasing n. 

For the p\up)ose of comparison with the exact solution, the equi- 
libriiim ten^ierature Tg for the flat plate can be obtained from 

equation (l)- 5 ) 'with B]_ = 0 (corresponding to q.w = O) 

^w = Tg = Constant. Thus, 


2o2 


(1 - 



2\t*1 


= 0 


or 



+ 


4^(1- 

a 



from which the recovery factor is 


Fr = 1 + 4 ^1 - a) 
a 
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The values of Fj, for six different tenrperature profiles are plotted 
against a In figure 2. These values are to he compared with the 

recovery factor from the exact solution ^taken as Fj- = which Is 
8l1so plotted In figure 2. 

The recovery factors resiiltlng from the use of the fifth-degree 
thermal profile (with 2 wall and 2 edge boundary conditions) are 
evidently acctirate to within a few percent over the range of Prandtl 
numbers applicable to gases • Note that for the same number of wall 
boundary conditions the recovery factor approaches unity as the nmber 
of stream conditions Is Increased. This somewhat anomalous result can 
probably be attributed to the effect of the assumption of equal velocity 
pnri thermal boundary-layer thicknesses. 

The wall shear stress 



Is Independent of the tenperature profile and can be written as 



by transforming to t]* and using equation (4o) for A and equations (l2) 
for ( — 3 :) with A = 0. 

wu 

Note that equation ( 50 ) can be put Into the same form as the 
corresponding resiilt In reference 12 by setting Cw = C \mder the 
Integreil only^ where C Is a constant equal to the average value of C^. 
This constant C was used In a similar manner In the solution of refer- 
ence 12; -hence, eqTiatlon ( 50 ) differs from the exact solution only In the 

value of the nxmierlcal factor V 57/515 which Is about 5 percent too hl^, 

Con 5 )arlson with stagnation flow; u^ approaches 0.- Approximate 
expressions for the heat transfer In the stagnation region of a blunt- 
nose body can be derived by noting that. In order to avoid Infinite 
values of d0/dx and dcp/dx as u^^ approaches 0 In equations (32) 
and ( 35 )^ "the following equations must be satisfied: 
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and 



P-wPw 12 + A 

(p\)^A 6 


9 


dx 




Bi = 0 


By using equation (15), these conditions can he written as 




= 0 


and 


£ - Tv = 0 

A ” oX 


(51) 


(52) 


Pt T T — 

where for this case — = = T^* For Un = 0, the quantities B2 

Pw -^1 

and B5 are both zero so that from equations ( 55 ) (56) 


^ = a + a^^Bq + ^(^ + o^Bl) 


A, 


( 53 ) 


and 



(5^) 


Note that the constants given in equation (^ 4 ) apply only for the partic- 
ular fifth-degree tenq)erature profile used in this section. For any other 
profile, the constants can he easily obtained from equation (25) and the 
desired boundary conditions. Using equations (55) and ( 54 ) and elimi- 
nating B]_ between equations (51) and (52) give the following relation 

between T^ and A: 



2k 
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12 + 
6A 




+ 


T. 


w 


2i 

A 


A 


2 

3 





( 55 ) 



( 22 ) and ( 23 ) are used for 


0 

A 


and 



Defining a Niosselt nmber for stagnation flow an 


Nu 




ki(Tv - t \) 


and vising eqimtion (31) gives 


Nu 



duj^ 

dx 



( 56 ) 


for constant Cp and a. This expression is plotted against T^/T^ in 

figure 3 for five tenperature profiles. The corresponding values from 
the exact solutions of references 13 and llf are plotted in figure 3 for 
conparison. Apparently, the present method can he expected to give 


T 

accurate results for this type of flow only for ^ and when the 

fourth- or fifth-degree temperatvire profiles are used with two edge con- 
ditions. Consldei*ation of the shape of the temperature and velocity pro- 
files given in reference I 3 indicates that the erroneous results of the 

present niethod for heating conditions [ — > 1 j are prohahly due to the 

\^1 / 

inherent limitations of the velocity profile rather than to the assunp- 
tion of equal thennal and velocity hovmdary-layer thicknesses. 


Comparison with experimental recovery factors: = 0.- The 

recovery factors on a flat plate were calcvaated with good accuracy 
by using the fifth-degree stagnation-temperature profile with two edge 
and two wall boundary conditions. (See fig. 2.) Therefore, as a 
further check on the range of applicability of the present method. 
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4 l > 


this particular profile only is used to con5)ute the recovery factors 
on a circular cylinder for comparison with the data of Eber (ref. 16). 

For zero heat transfer, q.^ = 0 ; hence, from equation (31) Bp = 0 
arirl the integral momentimi and energy equations (eq^s. (52) and ( 33 )) inay 
then he bitten as 



2 Z J2 + 


+ 




1 . ^ 
Up cLx 


2Cw 


^/Pl 0 12 + A 
Rl^p*lj A 6 


and 


(57) 


S t 5 ( 1 . 0 ( 58 ) 

dx N^l dx te - 1 dx / 


where 




t 


e 


and 




'Pi 


2 ^ 
T* 


tei-^ 

^1 dx 


( 59 ) 
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The quantity Z has "been introduced to avoid the use of. d^uq/dx^ as 
is done in the Holstein-Bohlen method. Equation (15) 'was used for 
in equation ( 59 )* 

A reiation between A and t0 can he obtained by integration of 
equation (58) which gives 

5 . ^ u-i (te - 1) = Consteint = 0 (60) 

since, at x = 0 , uq = 0 for a body with a blunt leading edge and 

^ = 0 for a body with a sharp leading edge. By using equations ( 35 )j 
Xj 

(5T)j (38), and (15), equation (60) can be written 


cr 

z 


12 + A 


dt. 


T*i q dUq ^ 

Tq u.q ^ 


a + ^ 

^ te2 + te 

tt 3 + I 


°-3 + I “7 


(1 - cr) 


7 “ 1 2 / A 

q2 + aP ^ 12 + A 


1 + 


[a^ + -^oj 



(61) 


which can be integrated by substituting a new dependent variable for 

te"^. Thus, equations (57), ( 59 ), (81) are to be solved simul- 

taneously for Z, A, and tg* 

In the terminology of reference 21 , the "local" recovery factor 
is defined as 


Fr = 



(62). 
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whereas a local free stream recovery factor is defined as 



(63) 


where in general Both of these recovery factors have heen 

calculated hy the present method for a circular cylinder with axis 
normal to the flow. The values of u^^ and du^^y^dx needed in the cal- 
culation were determined from impublished data obtained in the Gas Dynam 
ics Branch of the Langley Laboratory. These data were run at a free- 
stream Mach number of I.98 and free-stream Reynolds number based on cyl- 
inder diameter of 1.2 x 10^. The resulting pressure distribution could 
be closely approximated by the enrpirical relation 


£1 = 

Pco 


1VU^(0.67 + 0.98 cos 3-27 x) + 1 


(6l^) 


for 0 < ^ < 0.96. (The symbol L denotes the diameter.) 

The general procedxare \ised in the calculation was as follows: As 

a first approximation, it was assumed that te = 1; equations (57) and 
(59) were then solved simultaneously for Z and A by the method of 
isoclines . By using these values of A, a second approximation to the 
values of tg against x was obtained from equation (61). With the 
new values of te , the process can then be repeated to convergence . 
Actually, in the exa:iiq)le given, the second approximation for tg was 
sufficiently accurate since the use of this second approximation 
produced no significant change in the values of A against x 
determined from equations (57) and (59). 

The recovery factors Fj, and Fj,' resulting from this calculation 
have been plotted in figure Iv. The values of F37' from the data of 
Eber (ref. 16) and Fj. from the data of Eckert and Weise (ref. 22) are 
plotted in the .same fig\ire. Apparently, the conclusion alxeady obtained 
by several investigators for. subsonic flow (for exan^le, ref. 22 -id.th 

general discussions available in refs. 21 and 23) that F^ \fa, regard- 
less of body shape or axial station, is substantiated by the present 
results for supersonic flow about a cylinder. This result is also known 
to be valid for a pointed body of revolution as shown in reference 2 k. 
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Comparison vlth exact solution for heat transfer on a cylinder : 

T 

0 = 1; Ml «0: — .« 1 .-An exact solution for the heat transfer on a 

^ Tj_ 

cylinder was given hy Poland (ref. 15)- The solution postulates con- 
stant gas properties and hence woxild he strictly applicable only to 
incon^iressihle flow with T^/T]_ in the neighborhood of unity. As 
already phown by con^jeirison wlbh exact values (fig. 3 )^ the best restilts 

■T 

of the present method are in eiror by about 10 percent at ^ 1 for 

stagnation flow. Nevertheless, a con^iarlson of this method with Goland's 
exact solution would be useful as an indication of the accuracy of the 
present method with an arbitiary pressure gradient. 


For calculation of the heat transfer under the most general con- 
ditions, that is, in con^iresslble flow with arbitrary pressure and 
wall ten 5 )erature gradients and arbitrary Prandtl number near unity, 
it is convenient to write equations ( 32 ) and (33) in | the following form: 



1. ^ 
"^1 d 5 


2Cv 


R .2 

°L^1 \ 0 12 + 7v 

^1L\P*1/ ^ ^ 


( 65 ) 


and 


^ ^ + e fjL ^ ^ ^ /^1 L ?1 

dx I dSE \^1 (to Tv,- - 1/ oI^1l\Pi/ a h 


( 66 ) 


where 



(67) 
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and 



( 68 ) 


Equations (65) and (67) are of the same fonn as equations ( 5 T) ( 59 ) 

in the preceding section except that tg has heen replaced hy ^ and 
Bp 7^ 0 . Subtraction of the momentum equation ( 32 ) from equation (66) 
and use of equation (15) for a/L results in 


dx 


log H - 


5 * 


1 ) 


_1 
Mp dx 


^1 ^ 


- 1 



12 + a\ 


(69) 


which is more suitable for use in the general computational procedure 
than equation (66). The one-dimensional adiabatic energy equation 


du 


1 _ 


^1 


1 + 


7-1 




d% 


has heen used in equation (69) for the flow outside the boundary layer. 
Thus, equations (65), (67), (68), and (69) are to be solved simultaneously 
for Z, A, H, and Bp by a suitable modification of the method of 

isoclines . The quantities — , — , — , and 3 - are a1 1 functions of A 

A A A A 

and Bp as outlined previously. 

For application to the problem considered herein, equations (65), 
(69), (67)7 and (68) reduce to the form 


^ + 2 Z 
dx 


2 + y + (Tv - 1) 


10 - 5Bp 


1 


^ = 2(V 0 T2 + A 

i^l dx Rlj^ A 6 


(70) 
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log H 
dx 


1 + ^ + (Tt^ - 1 ) 


10 - 5Bi 

1 

-1] 

'"I 


du- 

dx 


1 _ 


T 1 ^"1 

1 i 

^ A 


12 + A' 


(71) 


z = 


^IL 



^1 



aoad 



a + 


®1®1 


|(ai^ + 


i^Bi) 




(72) 


(73) 


for B2 = B5 = 0, since = 0, a = 1, and Tyj- is a constant. Note 
that the gas properties are not necessarily constant. 

The same expression for the lix:^selt number as given by equation (56) 
can be used in this problem since T*^ = Tq_ = Tg = Too for low-speed 
flow. Thus, from equation (56) 


Nu _ ^ _ Mw / 1 ‘^^1 

koo(Tv - Tco) ypcoUooL 


(.1^) 


Inasmuch as the present method is expected to supply accurate heat- 
transfer values only under cooling conditions, the calculation has been 
T 

carried out for — = 0.7^- , This value gives starting values of A = 6.00 
and Bi = I.516 from equations (55) and. (52). The results in the form 
of NUy^^^ are shown, in figure 5 together with the corresponding values 
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given by Goland. Over the forward portion of the body the heat-transfer 
coefficients from the present method differ from the exact values by 
about 5 percent; however ^ at least half of this difference is due to 
the effect of variable gas properties as shown by comparison of Squire ' s 
solution (ref. lit) with the results of reference 13 • The assumption of 
equal thermal and velocity boundary- layer thiclcnesses apparently gives 
good results in this case. 

The confutation procedure for this particular problem was fairly 
sinfle since the energy and momentum equations were found to be practi- 

/ T 

cally independent . I For = 1 these equations would be entirely 

\ \ 

independent . ) Thus , as a first approximation^ B]_ was considered 

constant and equations (70) £Ln< 3 . (72) were then used to confute Z and 
X against x. With this curve of X against x, the quantities H 
and B2_ were confuted from equations (71) and ( 75 )- This second 
approximation for the Bn curve was then used to recompute Z and X 
from equations (70) (72); however, if allowance is made for confuta- 

tlonal error, the new values of X ane the same as the first values. 


APPLICATION TO FLOW WITH EQUILIBRIUM DISSOCIATION 


In a recent paper by Moore (ref. I7) the laminar boundary-layer 
characteristics with equilibrium dissociation were computed for a flat 
plate. One of the conclusions obtained was that the heat-transfer rate 
per unit area is essentially unaffected by dissociation when the plate 
tenferature is below 2 , 500 ° R. If an analogous conclusion could be 
obtained for the flow in the stagnation region of a blunt-nose body, 
the equilibrium dissociation would clearly be expected to have little 
effect on the heat transfer at any point (in the laminar flow) on an 
arbitrary body. Presumably, a heat-transfer calculation could be 
carried out on any body with known velocity and welU temperature distri- 
butions by the present method (as indicated in the section entitled 
General Method of Solution ). In view of the above statement, however, 
it is desirable to compute the heat transfer for the stagnation flow 
first. 

According to results previously given (fig. 5 ) Tor the heat transfer 
in the stagnation region with a perfect gas flow, there is little choice 
between the fourth- and fifth-degree temperature profiles for this pur- 
pose. Hence, for coii5)utatlonal sin5)llcity, a fourth-degree enthalpy 
profile with two edge and one wall boundary condition is used in this 
section. Thus, using equation (26) to express the density ratio Pq/p 
as a function of enthalpy and taking the limiting conditions from 
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equations (9) and (lO) for approacliing 0 (as indicated previously 

for a perfect gas) resiilt in the following two simultaneous equations 
for A and Bp: 


Kb + Kj^Bi + 


Pw 6A 


Pb 2 
Pq 315 




(75) 


and 



The K's are constants depending on the values of the stream enthalpy hQ_, 
the wall teniperature T^, and the h's in equation (26). 

Defining a wall Nusselt number in terms of the wall values and 
stream enthalpy as 



and using equation (31) for the heat transfer give a heat-transfer 
parameter of the form 


Nu^- 
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which is directly proportional to the heat-transfer rate for constant 
and h]_. This parameter Is plotted In figure 6 against h2_/hg for 

Hy = 2 , 000 ° R. The heat-transfer values for a perfect gas with no dis- 
sociation are also plotted In figure 6 together with the conrespondlng 
density changes . 

According to these resiilts, the equlllhrliim dissociation has a 
negUglhie effect on the heat- transfer rate at a stagnation point. 

This conclusion Is, of co\irse, to he regarded with caution owing to 
the Inherent l im itations of the method. Thus, for exan^jle, the velocity 
and enthalpy profile shapes are controlled only by stream nnfl -wall con- 
ditions . This limitation Indicates that application of the present method 
to the flat plate would result In no effect of dissociation as long as 
the stream and wall ten^jeratures are below dissociation values . Since 
this result agrees qualitatively with the results of Moore (ref. I7), 
the over-all profile shapes are apparently not too ing)ortant If the 
wall and stream values can be given accurately. 

Another disadvantage of the present method for application to dis- 
sociation flow Is associated with the assumption of equal thermal nnri 
velocity boundary-layer thicknesses. The dissociation data of refer- 
ence 17 show that the Prandtl number Increases to a TTiR-viTmim of 5 or 
6 for Increasing dissociation. This Increase would normally Imply a 
smaller thermal boimdary- layer thickness In cong)arlson to the velocity 
thickness. However, as noted previously, the con5)utatlons of reference 20 
Indicate that the thermal layer would be about 70 percent of the velocity 
layer even for an Increase In Prandtl number to 5 throughout the entire 
bou n da r y layer. Probably then, the present results are fairly reliable. 


CONCLUSIONS 


The application of the Von Karman-Pohlhausen Integral method to 
the laminar boundary layer of a flowing gas has been slD^jUfled for the 
general cane of arbitrary stream velocity and wall temperature distri- 
butions . The principal simplifying assumptions and procedures are 
(1) equal thermal, and velocity boundary- layer thicknesses, (2) linear 
viscoslty-teng)erat\ire relation, (5) tise of the first coefficient In the 
polynomial for the stagnation tenperature profile as one of the unknowns 
In the. final solution of the momentiim and energy equations, a.nr^ (Ij.) appli- 
cation of the Holsteln-Bohlen method to avoid the use of the second 
derivative of the stream velocity. 

Comparison of the resiilts of the present method with the Chapman- 
Rubesln nonlsothermal flat-plate solution (Jour. Aero. Scl., Sept. 194 - 9 ) 
shows that the sixth-degree (l wall boundary condition, 4- edge boundary 
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conditions), fifth-degree (2 wall 'boundary conditions, 2 edge houndary 
conditions), and sixth-degree (2 wall houndary conditions 5 edge 
houndary conditions) stagnation-ten5>erature profiles may he expected 
to yield accurate heat- transfer con^iutatlons except when the wall 
ten5)erature distrihutlons are characterized hy large greidients toward 
the rear of the hody, a condition which would seldom he realized in 
flight. The reason for the errors in the present method under this 
condition is prohahly the assun^ition of equal thermal and velocity 
boundary- layer thicknesses since the exact solution of Levy (Jour. 

Aero. Sci., May 1952) shows that the thickness of the thermal hoimdary 
layer is reduced hy increasing the wall ten5)erature gradient. 

On the other hand, the heat transfer in the stagnation region of 
a hlunt-nose hody is given accurately only hy the fourth-degree (l wall 
boundary condition, 2 edge boundary conditions) and the fifth-degree 
(2 wall boundary conditions, 2 edge boundary conditions) thennal pro- 
files when the wall teitperature is less than the local equilibrium 
ten5)erature . The error of the present method in this case, however, 
appears to he caused hy the inherent limitations on the velocity pro- 
file (as shown hy comparison with the results of Brown and Donoughe in 
MCA TN 214-79) rather than the ass4oii5)tion of equal thermal and velocity 
hound ary-layer thicknesses . 

The recovery factors on both the flat plate and circular cylinder 
in supersonic flow were con^juted satisfactorily with the fifth-degree 
(2 wall hoiindary conditions, 2 edge boundary conditions) profile. The 
local recovery factors on the cylinder were independent of location 
and nearly equal to the square root of the Prandtl number. These 
results are in agreement with those already obtained hy several inves- 
tigators for subsonic flow. 

The above res\ilts tend to indicate that the fifth-degree stagnation 
temperature or enthalpy profile will probably give the best over-all 
accuracy for the heat-transfer calculations by the present method on 
an arbitrary body for the case of heat transfer from the flowing gas 
to the surface. This conclusion is substantiated by conqparison of the 
results obtained from the use of the fifth-degree profile with the 
exact solution of Goland (Jour. Aero. Sci., July 1950) for "tlio heat 
transfer on a cylinder. 

As an illustration of possible applications of the method, the 
heat transfer was calculated in the stagnation region 'of a blunt body 
with equllibriTim dissociation. According . to this calculation the 
equillbri-um dissociation has a negligible effect on the heat-transfer 
rate at a stagnation point. However, the results obtained for this 
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extreme case should he taken as Indicative only of order of msignltude 
because of the inherent limitations of the method. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 19, 1953- 
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APPENDIX 

THE EFFECT OF LATERAL RADIUS OF CURVATURE ON SKIN FRICTION 
AND HEAT TRANSFER IN THE' LAMINAR BOUNDARY LAYER 


SelDEin and Bond (ref. 25) have given a solution by power series of 
the incon^ressihle -3 ami nar-houndary-layer equations for a slender 
circular cylinder with axis parallel to the flow. The conclusion was 
obtained that the changes in heat transfer and skin friction due to 
the effect of lateral ledius of curvature are negligible until the 
square root of the body Reynolds number is of the order of or less 
than 100 times the body fineness ratio. 

The results of this solution substantially confirm the earlier 
results of an approximate solution by Young (ref. 26) in regard to the 
changes in skin friction. Young used the Von Karm^-Pohlhausen Integral 
method but assumed the velocity-profile shape to be invariant over the 
entire body length. Althou^ the effects of lateral radius of curvature 
are certainly negligible for practical body configurations at subsonic 
speeds, these effects may become appreciable at supersonic speeds 
because of the large relative increases in wall tenqperature . 

The purpose of this appendix is to show how the application of 
the Von IttCriit^n-Pohlhausen method to this problem may be sin5)llfied 
for con^jressible flow by means of a modified Dorodnltsyn transforma- 
tion. The pressure gradient is assuned negligible since the results 
are likely to be of In^jortance only for bodies of high fineness ratio. 
However, in accordance with the physical requirements (even for zero 
pressure gradient) the dimensionless ten5)erature and velocity-profile 
shapes are allowed to change along the body length. 

The boundary-layer equations for the case of zero axial -pres s\ire 
gradient and constant Cp and, a are written as follows: 

Momentum: 


pur — + pvr — — = 
Sx oy 



(Al) 


0 = ^ 
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Continuity: 


Energy: 


^pur) + 


^pvr) = 0 

oy 




(A2) 


(A3) 


In terms of the stagnation temperature defined as 


Ql* 


T + 



eqiLiation (A3) may he written 



Integrating equations (A1) and (aJ 4-) across the boundary layer and 
using equation (A2) gives 


and 


(A5) 




1 *3Tv 
t*i dx 




dv = 

L ap^iU^LV^^- 


(a6) 
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If a modified Dorodniisyn transformation of the fonn 


% 





(A7) 


is introduced, equations (A5) and (a 6) can he written as 



and 



_ i _ r ^(1 

<3x Jq 




^'^/8t \ 


\ 'r/f^ 


(A9) 


The wall houndary conditions for the velocity and stagnation- 
teit5)erature profiles are derived from equations (A1) and (A^l-) evaluated 
at y = 0- When the viscosity and density relations given hy equa- 
tions (13) and (l4) and the transformation from equation (AT) are used, 
these wall conditions may he written as 


8^ \ ^2 Pw /8w \ 


(AlO) 


and 



2 _cos_p^/at \ 

P\\L J 


(1 


o) 


^1^ ftw 


(All) 
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Note that, for a = 1, these equations would he of the same form and, 
consequently, the velocity-profile shape w against T]/b and the 
stagnation-ten^erature-profile shape t against Tj/A would hoth be 
exactly the same if the same nijmber of edge boundary conditions were 
used in both profiles. Furthermore, if 6 = A (with a = l), equa- 
tion (A9) must reduce to the same form as equation (aS), so that under 

dT 

these conditions the additional equality — — = 0 is required. The effect 

dx 

of the -^mll temperature gradient for a body of high fineness ratio is 
probably negligible except on the forward part of the body where the effect 
of lateral radius of curvature is small. Thus, although the results are 
restricted considerably thereby, the use of the three assumptions /cr = 1, 

dT N V 

& = A, and — — = 0 ) should reveal the salient features of the effect of 

dx y 

lateral radi-us of curvature. 


Vnien these assuiig)tions are tised, the velocity and stagnation- 
teng)erature profiles are written 


t = w = 


(AI2) 


where, with three asyng)totic boiudary conditions at the edge^of the 
boundary layer. 


Ai = 


12 




\ 

\ 




\ 


+ A' 

/ aw ^ 

\ = 

/ a^ V 

6 

Wr. 

1 

\ r/w 

A' _ 1 


- 

“'' a2t “ 

2 2 


w 



w 




A5 = 


Ai, = 


3A* - 12 

6 

6 - A' 


(AI3) 




/ 
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The expression for A' is then obtained f ram either equation (AIO) or 
(All) as 


- 4^08 

cos Pvt A (\r 
3 L ■ ) 


(Al4) 


for 0 < a' ^ 12. 

Using equations (A12), (A13), and (A14) in 

gives 

^ ^ (12 + A* 

dx 18 uj^Lp*^ A' 


the momentum equation (a 8) 

(AI5) 


for a circular cylinder (r-w being constant). The same result would, of 
co\irse, be obtained from equation (A9) by using the assunq)tions discussed 
above. The value for 0 is obtained by integration of the expression 



Differentiating this equation and using equation (A14) for A gives 


5"^ - I 5TA' - - ^(A')5 , 

5 1A4 5 l44 

_ 12 + A' (12 + A' )^ _ 

(AI6) 

r^. Equating (AI5) and (A16) then gives a iinique rela- 
A ' and x as 


ix 315 l; 


for constant 
tion between 
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Pv^iL' 

57 A' 

5 1A4 



(12 + a *)5 


3T(V)^ - - 5(W 

3 i¥t- 

(12 + A')^ 


dA' 


Integration of equation (AI7) then results in 


(ait) 


5" = 


(-f 




A., 

61 

fV^lX 

35 

72 

Pw 



[37 a' 

1 

5 , 




23 


6(12 + A’ ) 


(12 + A')^ 


9 432' 


(A ')5 


+ f||log (12 + A’ ) - log 1^ > (A 18 ) 

(12 + A ’)^ ^1 


Thus, a' turns out to be a function only of the local body fineness 

, Pw^lX 

ratio x/ 2 rv and the wall Reynolds number ■ ■ ■ " . 

^^w 

The local 'skin- friction coefficient is defined as 


Cf = 


1 2 

|Pl^l 


or transforming to t)* and \ising equations (AI3) and {A±k) results in 


Cf = 4 




r 

iL 1 + A' Y 

I’w 6 / 


(AI9) 


The skLn-frictlon coefficient for the flat plate with constant wall 
temperature is 
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k2. 


and is derived from equation (50)> Thus, forming the ratio of the skin- 
friction coefficient of the cylinder to that of the flat plate gives 


(;, 20 ) 

^ \ 6 y »57 A'V6/ 


where A' is a function of | from equation (A 18 ) . Similarly, the 
local heat-transfer rate is 


Cf-r 


= 2 




or transforming to t]* and \ising equations (AI5) and (Al^t) give 




1 /rr,* \ 1 


18 


(A 21 ) 


"V 


The corresponding expression for the flat plate is 



since, with the ass\ui5)tions used herein, namely. 



the temperature profile on the plate is 


a = 1, 6 = A, and 


_T ^ % 

Ti Ti 


7-1 






2 
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(With the present assumptions, this ,same ten^jerature profile also satis- 
fies equations (A1) and (A3) for the slender cylinder.) The flat-plate 
relations from equations (l 2 ) and (i^O) with A = 0 and constant 

have also heen used to form . Forming the ratio of then 

P p 

results in 



(A22) 


The mean skin-friction coefficient is obtained hy integrating 
equation (AI 9 ) from x = 0 to x = L and using equation (AI 7 ) as the 
relation between dx and dA* . The result is 


Cf = ^ ^ 

105 L Pi 



12 + A* 
3 


^92 _ _ 5 _ (A*)^ 

12 + a' i1i4 12 + A' 


The mean skin-friction coefficient on the flat plate is 



The ratio of the mean skin-friction coefficient for the cylinder to that 
of the flat plate is 



1 /313 ^w 

210 V 37 L 



12 + A' 
3 


492 _ _ 5 _ (A')^ 

12 + A' ll<4 12 + a’ 


(A23) 


where A' must correspond to the values from equation (A18) with x = L. 
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In figure J, equations (A20) and (A25) are plotted against 

which is similar to the parameter used by Sehein and Bond- 






(ref. 25). The ratio from reference 25 is plotted against 


I' 


in figure 7. 


Figure 7 shows that the results from the present method are in 
good agreement with those from reference 25 if the stream Reynolds 
number ppuj_x/p.3_ in reference 25 is replaced by the wall Reynolds 

nimiber The con^ressihility effect on the slender cylinder 

is thus approximately accounted for hy evaluating the gas properties 
at the wall rather than in the free stream. This result might have 
been anticipated in view of a discussion in reference 27 on the com- 
pressibility effect in the 1 aminar boundary layer. 


The effects of the lateral radius of curvature on both the skin 
friction and heat transfer are increased over the incoii5)resBible effects 
in almost direct proportion to T^/T3_ since 


Tl 


with the present assumption for the viscosity variation. As a practical 
exan^jle, consider a body of fineness ratio I5 at Mp = 5*00 and 

RIt = 2 X 10^ irLth — = 5* Under these conditions the value of 
at the rear of the body is 


II 


L/^w 



0.11 


(|j^ is the value of | where x = L.) Tben^ from figure 7 the local 

sld.n-friction and heat- transfer coefficients would be about 20 percent 
greater than those on a flat plate for the same stream conditions. The 
meein slcLn-frlction coefficient and, incidentally, the total heat transfer 
would be about 7 percent gi*eater than on the plate. 
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TAET.B I.- NUMERICAL VALUES OF a FOR EQUATION (55) 


Thermal profile 

a 

ai 

02 

05 


05 

06 


Degree of 
polynomial 

Number of 

boundary 

conditions 

Wall 

Ed^e 

5 

2 

2 

0.581746 

- 0.146052 

-0.055175 

-0.013889 

0.044841 

-0.012302 

-0.003572 

-0.000794 

6 

2 

5 

.29805 

-.09582 

-.028067 

-.00552 

.041254 

- .010157 

-.002489 

-.000453 


0 

3 

.17381 

- .028175 

0 

0 

.029762 

-.005968 

0 

0 


1 

2 

.326190 

-.104565 

-•025597 

0 

.041667 

-.009921 

-.001984 

0 

5 

1 

5 

.24286 

-.062698 

-.011508 

0 

.036905 

-.007540 

-.001190 

0 


1 

4 

.187662 

- .040625 

-.005988 

0 

.032576 

-.005808 

-.OOO75Q 

0 
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dT 

Relative heat -transfer values with =* 0; n =• 0 

Degree of 
thermal 
profile 

Number of 
boundary conditions 

1 Description of thermal profile 

Yn’(O) 

— 

Percent 

error 


Edge 

Degree of 
polynomial 

Number of 
conditions ' 
at wall 

Number of 
conditions 
at edge 

O 4 
D 5 

C> 6 

□ 4 
V 6 
> 5 
A 4 

H 

1 

2 

3 

3 

4 
3 
2 

4 

4 

5 

6 

4 

5 

6 

0 

1 

1 

1 

2 

2 

2 

3 

2 

3 

4 
1 
2 
3 

0.543 

.601 

.577 

.526 

.541 

.575 

.576 

-8.2 

1.6 

-2.5 

-11.1 

-8.5 

-2.8 

-2.6 

\ 

Exact, ref. 12 

0.5915 

0 



Figiire 1 .- Variation of heat-transfer parameter with 
isothermal flat plate with cr = O.72. The symhol 







n 

n 


for the non- 
is defined hy 


the equation 










































Recovery factors, 
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A 2.50 8.68 X 1Q4 

□ 1.86 1.49 XlO^ 

O 1.56 1.88 X 10^ 


Data of 
reference 16 


Present method -using 

equation (64) for 

■with M„ = 1.98, CT = 0.72 


-Fj. From data of 
reference 22 


Distance from stagnation point, x/L 

Figure 4.- Local free-stream recovery factor Fj.' and local recovery 
factor Fr on a circular cylinder -with axis normal to the flow. 
Separation occurs at approximately — = 1.0 where L Is the 


Is the 


cylinder diameter. 






Heat-trcinsfer 



Figure 5.- Heat transfer on a cylinder the present method and from 

Golan d ' s exact solution. 


U1 

VjJ 
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No HlRRon.lfltlnn T.„ - PnOnO R.- rr = 0 flR 

J - yy - - J 

Dissociation at 0.001 atmosphere p 
T^- 2000°R, o^ = 0.65 


Heat-transier parameter 


lO wSjj. ueuaity ra.LiO 


Local stream enthalpy parameter, 


Figure 6.- Variation of heat- transfer parameter with stream- enthalpy- 
parameter with and without diBsoclation at the stagnation point of a 
blunt body. Wall temgperature, Ty = 2,000° R. 
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